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ABSTRACT 


The  Sexico  v^o  resistance  and  self -propulsion  model  test  results  are 
represented  by  curve  fitting  methods  which  result  *n  functions  containing 
Chebyshev  polynomials.   The  accuracy  of  the  curve  fitting  method  is  dis- 
cussed, and  the  range  of  validity  of  the  functions  is  presented. 

A  computer  program  for  the  calculation  of  EHP  and  SEP  for  single- screw 
merchant  ships  with  hull  proportions  falling  within  the  range  of  the 
Series  60  is  discussed  in  detail.   Several  examples  in  the  use  of  the  pro- 
gram are  presented. 

It  is  concluded  that  estimates  of  EHP  and  SHP  can  be  calculated  to 
within  a  tolerance  in  the  neighborhood  of  1$  on  the  average,  and  with  great 
certainty,  to  within  a  tolerance  of  less  than  Zl°  for  any  given  speed. 
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I.      INTRODUCTION 

1.  Foreword 

The  determination  of  the  required  horsepover  necessary  to  fulfill 
the  speed  requirements  of  a  new  ship  design  is  one  of  the  many  problems 
vhich  faces  the  naval  architect  in  the  early  stages  of  ship  design. 
He  must  have  some  way  of  knowing  what  effect  various  choices  of  size, 
proportions,  and  fullness  of  the  ship  -will  have  on  the  power  required  to 
propel  the  ship  at  the  speeds  he  desires. 

In  determining  this  pover  requirement,  the  naval  architect  can  drav 
upon  his  past  experience  in  ship  design  and  the  results  of  specific  isolated 
model  tests  which  are  published  periodically  in  the  various  technical 
journals  of  the  profession. 

However,  thr   results  of  experiments  conducted  on  families  of  models, 
in  which  the  different  design  parameters  are  varied  systematically,  are  the 
naval  architects'  most  valuable  sources  of  reference  in  estimating  the 
power  required  for  his  new  design.   The  result  of  these  tests  are  usually 
expressed  in  the  form  of  design  charts  from  which  the  naval  architect, 
using  interpolation  techniques  where  necessary,  can  calculate  the  power 
requirements  for  a  number  of  combinations  of  design  parameters,  and  thereby 
make  an  intelligent  choice  of  the  best  combination  of  parameters  to  give 
minimum  power  with  due  consideration  for  the  other  limitations  of  his  design 
conditions.   The  Series  60  is  such  a  family  of  models  designed  to  aid  the 
naval  architect  in  the  estimation  of  power  requirements  for  single-screw 
merchant  ships. 

2.  The  Series  6o  Model  Tests 

The  research  on  Series  60  was  carried  out  at  the  David  Taylor  Model 
Basin  (DTMB)  of  the  United  States  Navy.   Experiments  were  conducted  as  part 
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of  the  Bureau  of  Ships  Fundamental  Hydromechanics  Research  Program  during 
the  years  19*+8  to  i960.   The  results  were  published  in  a  series  of  papers 
presented  before  the  Society  of  Naval  Architects  and  Marine  Engineers, 
and  a  final  report,  incorporating  all  of  the  information  contained  in  the 
interim  reports  was  published  in  July  of  1963.   ; 

After  a  survey  of  current  practice  in  the  shipbuilding  industry  to 
ensure,  as  far  as  possible,  that  the  series  would  cover  the  normal  range 
of  proportions  of  modern  ships,  proportions  were  chosen  for  five  parent 
models  of  the  Series  60.   These  proportions  are  listed  in  Table  I  below. 

TABLE  I 
Hull  proportions  for  Five  Parent  Models  of  Series  60 


Model  No. 

U210 

1+218 

-1+221 

U213 

1+2  Ik 

cb 

.60 

.65 

•  70 

•75 

.80 

l/b 

7.50 

7.25 

7.00 

6.75 

6.50 

b/h 

2.50 

2.50 

2.50 

2.50 

2.50 

A/(l/ioo)3 

122.0 

lto-u 

I63.I+ 

188.2 

216.5 

L/71/3 

6.I65 

5.869 

5-593 

5.335 

5.092 

LCB  as  fo  of  L 

from  © 

**  BP 

1.5  aft 

I.5I+  aft 

0.55  aft 

1.5  fwd 

2.5  fwd 

For  any  one  block  coefficient  and  LCB  position,  a  total  of  nine  models  were 
tested  in  which  the  L/B  and  B/H  ratios  were  varied,  for  a  total  of  forty-five 
models.   To  cover  the  general  spread  of  L/B,  B/H,  and  A/(L/l00)  for 


Numbered  references  throughout  the  text  are  listed  on  rage  h  h 
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existing  designs  and  the  possible  variation  in  LCB  position,  a  grid  was 
adopted  as  shown  by  the  dotted  lines  in  Figure  I. 

* 

FIGURE  I 

Variation  of  Proportions  with  C^  for  Series  60  Models 


*  Frci^  reference  (l),  page  II- 3 
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The  pattern  for  a  typical  case  (C,  =  .60)  is  shown  in  Figure  II, 

FIGURE  II  * 

Typical  Variation  of   L/B  and  b/H  Ratios  for 
a  Given  Value   of   C, 


3.5 


© © 


65 


|x  3.0  (t) (t) 


s© 6 


7.5 

e 


© 


© 


8.5 


The  midship  section  of  the  node Is  had  no  deadrise  in  accordance  with 

current  practice.   The  relationship  "between  C  and  C  for  the  Series  60  is 

m      p 

as  follows: 

c  =  0.0657  c  +  .925 

m  p 


And,  since    C,  =     C  x  C 
b        p    m 


CL  =  O.O857  C   +  .925  C 
b  P         p 


(1) 


LWL  for  the  Series  60  hull  forms  is  related  to  L  by  the  following 

simple  formula : 

LWL  =  UO0.75  x  L  (2) 

400. 

*  From  reference  (l),  page  II-3 
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The  results  of  the  Series  60  resistance  tests  are  presented  in  the 
form  of  easy  to  use  design  charts,  and  in  numerical  point -by-point  form 
in  tables. . 

Values  of  total  resistive  coefficient  (C, )  for  forty-five  ships  of 

length  betveen  perpendiculars  of  k-CO  feet,  and  varying  values  of  C  ,  L/B, 

b 

and  E/H  are  presented  in  Tables  B-l  through  B-l+5  of  reference  (l). 

Utilizing  the  Froude  assumption  that  the  total  resistance  can  be  divided 

into  two  parts,  the  skin  friction  of  an  "equivalent  plank"  and  the 

residuary  resistance,  the  values  of  C,  include  a  skin  friction  coefficient 

(C  )  calculated  for  salt  water  at  59°?  in  accordance  with  the  American 

Towing  Tank  Conference  (ATTC)  19^-7  line  ,  and  an  additional  allowance  for 

ship  oo^relation  (AC  )  of  +  0.000^  as  recommended  by  the  ATTC  I9U7. 

Methods  for  scaling  the  values  of  C  for  an  arbitrary  ship  length  will  be 

discussed  in  Chapter  TV. 

Values  of  propulsive  coefficient  (PC  or  EHP/SE?)  for  forty-five  ships 

of  length  betveen  perpendiculars  of  6C0  feet,  and  varying  values  of  C  , 

b 

L/B,  and  B/n  are  also  presented  in  Tables  B-l  through  B-l+5  of  reference  (l). 
Unlike  the  values  of  C  ,  the  propulsion  data  can  not  be  corrected  for 
variation  in  ship  length.   The  choice  of  a  600  foot  ship  to  illustrate  the 
propulsion  tests  vas  made  principally  because  it  vas  considered  more  rep- 
resentative of  modern  ships  than  the  h-OQ   foot  length  chosen  for  the  resistance 
presentation.    The  propulsion  data  presented  vas  obtained  using  a  standard 
propeller  diameter  equal  to  O.70  of  the  designed  lead  draft. 

Data  for  the  effect  of  variation  of  LCB  position  on  power  requirements 
is  by  no  means  complete.   Corrections  to  total  resistance  for  variation,  of 


* 

Also  known  as  the  Schcenher  line. 

Reference  (l),  page  VII- h. 


-5- 


LCB  position  for  five  values  of  C  are  presented  in  Tables  ky   through. 

b 

53  in  reference  (l),  with  only  four  or  five  different  positions  of  LCB 
noted. 

Although  the  propulsive  data  is  presented  for  ships  of  600  foot 
length,  if  the  designer  adopts  the  lines  of  Series  60,  the  position  of 
LCB  as  used  in  the  parent  forms,  and  a  propeller  having  the  standard 
ratio  of  diameter  to  draft  of  0.7,  he  can  make  a  very  accurate  estimate  of 
both  the  EHP  and  3H?  of  a  ship  for  any  particular  selection  of  L,  B,  H, 
and  A  . 
3.   Previous  Work  in  Series  60  Computerization 

The  development  of  a  computer  program  for  estiniating  EHP  using  the 

(2) 
Series  60  results  was  first  attempted  by  D.  Khoushy  in  1962. *    This 

attempt  represented  Lhs  Series  60  by  curve  fitting  methods  to  some  extent 

e.g.  wetted  surface,  but  it  relied  exclusively  on  interpolative  methods 

for  C  which  are  both  cumbersome  to  program  and  highly  consumptive  of 

at 

computer  storage  space  .  Furthermore,  the  large  amount  of  storage  space 
required  in  this  type  of  program  made  it  unsuitable  for  use  as  an  EHP  or 
SHP  subroutine  in  a  larger  program.   It  also  should  be  noted  that  an 
interpolative  scheme  is  not  very  efficient  since  the  test  results  are 
not  presented  in  a  uniform  manner,  i.e.  at  low  values  of  C,  models  were 


tested  up  to  a  V/VlWL  of  1.10  whereas  the  higher  values  of  C  were  tested 
up  to  a  V/^fL"WL  of  0.80.   This  problem  was  avoided  in  the  Khoushy  work  by 
discarding  a  large  amount  of  Series  60  data  in  the  low  and  high  ranges  of 

v/YuTl. 


The  storage  of  approximately  6C00  numbers  is  required  for  the  Series  60. 
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Therefore,  it  was  decided  that  curve  fitting  vculd  be  the  most 
advantageous  method  of  representing  the  Series  60  model  test  results 
in  a  computer  program. 
k.      The  Curve  Fitting  Approach 

The  curve  fitting  approach  for  representing  the  Series  60  results 
is  very  desirable  since  once  a  function  is  found,  which  represents  the 
given  data,  one  knovs  the  error  he  can  expect  in  predicting  EHP  and  SHP 
using  these  derived  functions.  The  fitting  of  curves  through  the  Series 
60  resistance  data  was  started  in  early  1966  by  T.  M.  Pitidis-Poutous. 
The  basic  theory  of  the  curve  fitting  method  employed  was  developed  in 
reference  (3)  in  19&5  "with  the  mathematical  derivation  of  ships  lines 
being  the  focus  of  attention.   This  curve  fitting  method,  as  it  applies 
to  the  representation  of  the  Series  60  resistance  data  and  its  extension 
to  represent  the  propulsion  results  for  this  work,  is  described  in 
Chanter  II. 
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II.   DESCRIPTION  OF  CURVE  FITTING  METHODS  USED  FCR  THB  SERIES  60 

The  theory  of  the  curve  fitting  methods  used  for  representing  the 
Series  60  model  test  results  stated  herein  is  largely  adopted  from  reference 

w. 

Since  the  chart  reading  methods  of  the  Series  60  are  not  suitable  for 
computerization,  it  vas  decided  to  fit  curves  through  the  numerical  results 
presented  for  C  and  PC  (or  EHP/SH?)  in  Tables  3-1  through  B-U5  of  reference 

(1). 

In  the  Series  oO,  both  C     and  PC  are  functions  cf  four  independent 


variables,   i.e.   V/YiTTl,    C  ,    L/B,   and  B/H.     Thus,   either  C     or  PC  can  be 

b  t 

represented  by  an  equation  of  the  form: 


y=|  tit  4nhji  ^)%(^)xkK)aJ^    0) 

and  the  problem  reduces  to  determining  the  matrix  of  coefficients  1^-wkJ.i.J 
from  the  given  data.   If  the  coordinate  functions  Q  ,  ID  ,  }(,  and £1  are 
taken  as  ordinary  polynomials,  the  above  expression  becomes 

a«i  j-i  k=«  n=i 

It  can  be  shown  that  there  are  certain  advantages  to  be  gained  if  the 

coordinate  functions  are  taken  as  the  Chebyshev  polynomials,  with  no 

additional  complications.    The  independent  variables  V/TEWL,  C  ,  L/3, 

b 

See  reference  (5). 
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and  B/H  can  be  normalized  in  the  range  of  -1  to  +1  by  the  equations  given 

in  Table  II,  and  the  corresponding  symbols  used  for  the  original,  normalized, 

and  computer  program  variables  are  given  in  Table  III.  The  equations  of 

Table  II  vere  derived  on  the  basis  of  the  limits  of  the  Series  60  model 

tests. 

Stipulations  on  the  maximum  powers  I,  J,   Kf   and  II,  to  be  used  in 

equation  (3)  are  determined  by  the  amount  of  data  available  from  the 

Series  60  results.   For  example,  each  model  was  tested  for  three  different 

values  of  B/H  and  L/B,  therefore  I  and  J  must  be  equal  to  or  less  than  3* 

Models  were  tested  for  five  different  values  of  C  ,  therefore  K  must  be 

b 

equal  to  or  less  than  5*   Thus,  the  number  of  terms  in  the  equation  will 
be  equal  to  or  less  than  ^5N.   The  value  of  N  cannot  be  arrived  at  in  the 
same  way  as  I,  J,  and  K,  since  each  model  was  tested  for  twenty  to  thirty 
different  values  of  V/TfLWL,  differing  from  model  to  model.   If  we  assume, 
for  the  moment,  that  N  =  10  would  be  adequate,  we  see  that  the  determina- 
tion of  the  coefficients  OL^^-   £   would  involve  the  solution  of  ^50 
linear  algebraic  equations.   This  is  virtually  impossible,  both  from  the 
numerical  and  computer  storage  points  of  view. 

These  considerations  prohibited  the  use  of  a  direct,  if- parameter  fit, 
and  indicated  the  necessity  of  a  slep-by-step  technique. 

Using  the  symbols  for  the  normalized  variables  as  defined  in  Table  III, 
equation  (3)  becomes: 

or   *   •» 


y*ZT.  A«*  5 (*■**)  nn(*) 


(*) 


i-i  jt-i 
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table  ii 


Equations  for  ITorrali zing  Independent 
Variables  in  Series  50  Ran^e 


v  =  2.5    _v—      -  1.875 

YTJE 


-L-  =  O.kv  +  0.75 

Vlwl 


u  =     L/3       -  7.0 
1.5 


E 


=  1.5u  +  7-0 


=  2.0     B        -  6.0 


H 


B 


H 


o.5z  +  3.0 


x  =  10.0     C        -   7.0 


=  O.lx  +  0.7 
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vhere 


A  =  j    +    3"(x-/)    +    IJ  Ck-0 


and  F,     -        $;(*)      fyly)     Xk(*)        . 

The  above  substitutions  reduce  the  f cur-dimensional  matrix  "a"  to 

a  two-dimensional  matrix  "A",  which  is  easier  for  computations,  and 

defines  the  problem  more  appropriately  for  the  method  to  be  used.   For 

given  values  of  u,  z,  and  x,  equation  (h)   can  be  written  as: 

N 


yp  -  £  c„  n,  h 


(5)  . 

with  L. 

Qn  »  Y.      A^  ^  ("'>'^  (6) 

-*  =  " 

Equation  (5)  represents  the  variation  of  Y  with  respect  to  v  for 
the  pth  model,  and  the  N  coefficients  (Cpa,  n  =  1,  N)  can  be  evaluated 
from  the  data  in  Tables  B-l  through  B-^5  of  reference  (l),  using  the 
least-squares  criterion: 


v/n 


by  differentiating  partially  with  respect  to  the  unknown  coefficients. 

All  models  were  fitted  in  this  fashion,  and  the  resulting  equation 
for  each  model  was  evaluated  at  intermediate  values  in  order  to  ensure 
that  the  polynomial  had  no  undesired  oscillations.   In  fitting  the  curves 
for  C  ,  it  was  found  that  a  reasonably  good  fit  was  obtained  with  N  =  10 
(corresponding  to  a  9th  order  polynomial)  and  since  this  is  rather  high, 
the  possibility  of  introducing  a  suitable  weighting  function  was  explored. 
Therefore,  the  fittings  were  repeated  with  the  following  two  alternative 
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weighting  functions  assigned  to  C  : 


Y=    (Ct*i*')(>//YZZL)S 


The  last  alternative  was  found  to  give  the  best  fit  with  N  =  10  for  most 
models.   Going  from  N  =  10  to.  N  =  11,  the  maximum  error  experienced  in 
C  x  10^  only  decreased  from  .025  to  .021+  for  Table  B-l  of  reference  (l). 

Having  obtained  matrix  £c  "\  corresponding  to  equation  (5)  with 
p  =  "1+5  (number  of  models),  matrix  fA  ~\   was  obtained  by  exact  fit,  i.e. 
L  =  1+5*   It  should  be  mentioned  here  that  a  certain  amount  of  regression 
could  also  be  tried  on  this  step  of  fitting,  and  this  was  actually 
attempted,  but  without  success.   The  reason  is  that,  although  it  might 
be  possible  to  obtain  a  good  fit  for  C  or  PC  at  a  constant  value  of 
V/VlWL  with  L  less  than  1+5,  it  does  not  follow  that  such  a  fitting  could 
be  applied  successfully  to  the  coefficients  of  the  C  -  v/VLWL  equation 
(i.e.  C   ).   Thus,  direct  U-parameter  fitting  of  the  data  could  have 
produced  an  equation  with  less  error  and  less  terms,  if  it  were  possible 
numerically. 

The  values  of  PC  were  fitted  in  the  same  manner  described  above,  but 
no  weighting  function  was  used. 

The  final  equations  obtained  are  as  follows : 

,!£=-  'tl     A-^K^T.M  (7) 


10       4T 


PC=KMX.  b^k>^t„.,m 


(8) 
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where 


I      =  j  +  3(x-i)  +  ?(*-i) 

jp   =.  p  *k     Chebyshev  Polynomial 

The  elements  of  the  ma  trice-  \  A  land  |B  "\  are  read  in  as  data 

1  nl1    l  nlJ 

in  the  main  program,  and  are  listed  in  Appendix  A.   These  matrices  form  a 
compact  computational  package  of  900  numbers  for  the  calculation  of  C 
(400  ft.  ship)  and  PC  (600  ft.  ship). 

No  attempt  was  made  at  fitting  curves  through  the  data  for  correction 
to  total  resistance  due  to  variation  of  LCB  position.   The  inadequacy 
of  the  data  presented,  as  veil  as  the  erratic  nature  in  which  it  behaves 
with  v/Yu.rL,  precluded  the  use  of  the  curve  fitting  methods  described  above 
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III.  RESULTS  OF  CURVE  FITTING  METHODS  USED  AND  RANGE  0?  VALIDITY 

The  results  of  the  curve  fitting  methods  used  for  C  and  PC  were 
first  tested  by  reproducing  the  date  contained  in  Tables  B-l  through  B-^5> 
i.e.  for  the  forty-five  models  of  the  Series  60.   The  results  of  this 
test  are  presented  in  Table  IV,  where  the  maximum  error,  in  percent, 
experienced  in  C  and  PC  is  noted  for  each  of  the  forty-five  models. 

The  maximum  error  experienced  in  C  was  2.7^0$  for  the  parent  model 
1*2 18.  The  maximum  error  experienced  in  PC  was  2.030$  for  model  ^276. 
Although  these  figures  appear  prohibitive  at  first  glance,  it  should  be 
emphasized  that  these  are  the  maximum  errors  experienced  for  any  given 
model,  and  they  occur  at  only  one  value  of  v/{jLWL.   The  errors  at  other 
values  of  v/VLWL  for  any  given  model  are  well  under  Vfc     in  the  great 
majority  of  cases.   In  fact,  as  can  be  seen  from  Table  IV,  the  maximum 
errors  experienced  in  C  and  PC  for  any  given  model  are  well  below  Vfo 
in  most  cases.   To  illustrate  this  point  further,  plots  of  C  versus 


V/YLT>7L  are  presented  in  Figures  III  through  VII,  and  plots  of  PC 
(or  EH?/SEP)  versus  V*/VlWL  are  presented  in  Figures  VIII  through  XII, 
for  the  five  parent  models  of  the  Series  60.   In  these  plots,  the  solid 
line  represents  the  value  of  C  and  PC  as  calculated  from  equations  (7) 


and  (8)  respectively,  and  the  points  represent  the  value  at  that  V/VLWL 
presented  in  the  tables  of  reference  (l),  which  are  in  error  enough  to 
be  visually  discernable  on  the  scale  used  in  making  the  plots.   It  is 
further  noted  that  Figure  IV  represents  the  worst  case  experienced  in 
fitting  either  C  or  PC. 

As  a  second  test,  equations  (7)  and  (8)  were  evaluated  by  generating 
125  interpolate  and  extrapolate  models,  especially  concentrating  on 
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TABLE 

IV 

Maximum  ' 

)c  Errors  Experienced 

in  Ct  and  PC  for 

the  45 

Model?   of  the 

Series  60 

Series 

60       Model 

] 

)tex.  dp  Error 

Max  $  Error 

Taole  ] 

No .            No . 

h2ko 

.60 

l/b     b/h 

6.5    2.5 

in  ct 

.578 

in  PC 

B-l 

•853 

B-2* 

4210 

.60 

7-5      2.-5 

1.230 

.648 

B-3 

4267 

.65 

8.25     2.5 

1.515 

.672 

B-4 

4264 

.65 

6.25     2.5 

.340 

.616 

B-5* 

4213 

.65 

7.25     2.5 

2.740 

.298 

B-6 

4243 

.60 

8.5       2.5 

.701 

.696 

B-7 

4244 

•  70 

6.0      2.5 

•  531 

.  .248 

B-8* 

4221 

.70 

7.0      2.5 

.370 

.402 

B-9 

4247 

.70 

8.0      2.5 

.232 

.552 

B-10 

4268 

•  75 

5-75    2.5 

.819 

.336 

B-ll* 

4213 

•  75 

6.75    2.5 

.824 

.697 

B-12 

4271 

•  75 

7-75    2.5 

1.045 

.682 

B-13 

4243 

.80 

5-5      2.5 

.191 

.372 

B-l4* 

4214 

.80 

6.5      2.5 

.331 

.134 

B-I5 

Wdirj 

n  q      2.5 

.243 

.394 

B-l6 

4252 

.60 

6.5      3.0 

.885 

.437 

B-17 

.  4255 

.60 

7.5      3.0 

.714 

•  393 

B-13 

4254 

.60 

8.5      3-0 

.477 

I.261 

B-19 

4272 

.65 

6.25    3.0 

.602 

.529 

B-20 

4275 

.65 

7.25    3-0 

.405 

.446 

B-21 

4274 

.65 

8.25    3.0 

.600 

.516 

B-22 

4256 

.70 

6.0      3-0 

.175 

•     .277 

B-23 

4259 

.70 

7.0      3.0 

.134 

•  590 

B-24 

4253 

.70 

8.0      3-0 

.278 

•  573 

B-25 

4276 

.75 

5-75     3.0 

.528 

2.030 

B-26 

4279 

.75 

6.75    3-0 

1.331 

1.430 

B-27 

4273 

•75 

7.75    3-0 

.923 

.591 

B-23 

4260 

.80 

5-5      3-0 

.243 

.517 

B-29 

4263 

.80 

6.5      3.0 

.166 

.129 

B-30 

4262 

.60 

7-5      3-0 

.082 

.284 

B-31 

4241 

.60 

6.5      3.5 

.5^1 

.449 

B-32 

4253 

.60 

7-5      3.5 

.411 

.878 

B-33 

4242 

.60 

3.5      3.5 

.459 

.906 

B-34 

4265 

.65 

6.25    3-5 

.634 

•  537 

B-35 

4273 

.65 

7.25    3-5 

.418 

.349 

B-36 

4266 

.65 

8.25    3-5 

.493 

.664 

B-37 

4245 

•  70 

6.0      3-5 

.203 

.443 

B-33 

4257 

•  70 

7.0      3-5 

.435 

.452 

B-39 

4246 

.70 

8.0      3.5 

.312 

.544 

B-4o 

4269 

•  75 

5.75    3.5 

..223 

.598 

B-4l 

4277 

•  75 

6.75    3-5 

.506 

.581 

B-42 

4270 

•  75 

•7-75    3-5 

.919 

1.275 

B-43 

4249 

.80 

5.5      3.5 

.192 

.542 

B-44 

4261 

.80 

6.5      3.5 

.191 

•  577 

B-45 

4250 

.80 

7.5      3-5 

.118 

.303 

*  Denotes  parent  model 
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intermediate  values  of  the  fiveCb's  used  for  the  Series  60  models.  This 
vas  done  in  order  to  determine  the  range  of  validity  of  these  equations 
with  respect  to  their  four  independent  variables.   It  should  be  emphasized 
that  the  ranges  of  validity  determined  below  are  not  just  a  result  of 
the  inadequacies  of  the  curve  fitting  methods  used;  they  are  very  much  a 
result  of  the  limits  of  the  Series  60  model  tests  themselves  and  the 
range  of  data  made  available  in  reference  (l). 

The  data  presented  in  reference  (l)  is  based  on  the  use  of  three  values 
of  B/H,  i.e.  2.5,  3*0,  and  3«5«   Since  equations  (7)  and  (8)  rely  on  a 
second  order  interpolation  scheme  in  B/H,  which  is  the  method  suggested  in 
reference  (l)  when  using  the  design  charts,  there  should  be  no  discrepancies 
caused  by  this  variable.   It  is  also  expected  that  a  limited  amount  of 
extrapolation  is  possible  in  B/H. 

Similarly,  there  should  be  no  discrepancies  due  to  L/B  since  equations 

(7)  and  (8)  rely  on  a  second  order  interpolation  scheme  for  this  variable^ 

which  is  also  the  method  used  by  BTI-B  in  the  preparation  of  the  Series  o0 

design  charts.  However,  as  can  be  readily  seen  from  Figure  I,  there  are 

limits  (imposed  by  the  proportions  of  the  Series  60  models)  on  the  range 

of  interpolation  possible  in  L/B.  The  range  of  interpolation  in  L/3  is  a 

function  of  C  and  can  be  expressed  as  follows: 
b 


(L/B)^  =  5-5+    8-°  -  1Q  Cb  (9) 

2 


'min 


(L/3)^  =   7.5  +    8.0  -  10  Cb  (10) 

2 
_ 

See  example  of  the  ship  SCHUYLER  OTIS  BLAirD  in  Chapter  V. 


-25- 


A  small  amount  of  extrapolation  "beyond  these  limits  should  not  cause 

trouble . 

In  producing  the  Series  60  design  charts,  a  third  order  interpolation 

in  C  vas  used  by  DUMB,   In  equations  (7)  and  (8)  a  fourth  order  inter- 

polative  scheme  for  C  is  used  vhich  could  produce  slight  discrepancies 

b 

in  interpolation  and  extrapolation.   It  is  recommended  that  the  program 

"be  used  only  in  the  range  of  C  =  0.60  to  0.80. 

b 

The  Series  cO  model  tests  impose  definite  limits  on  the  range  of 

V/VLWL  permissible.   For  example,  in  testing  equations  (7)  and  (8)  with 

the  125  models  mentioned  above,  it  vas  found  that  with  C  =  O.775,  b/h  =3.0, 

and  L/B  =  7.0,  the  values  of  C.  and  PC  given  by  equations  (7)  and  (8) 

dropped  off  sharply  at  V"/VlWL  =  .925  and  even  went  negative  at  v/YUfl  =   .950. 

This  is  understandable  since  at  these  values  of  v/TLWL  ve  are  outside  the 

range  of  results  given  by  the  Series  60  tests.   Therefore,  to  avoid 

occurrences  of  this  nature,  a  "trap"  has  been  built  into  the  EHP-SKP  program 

described  in  Chapter  IV,  vhich  puts  limits  on  the  value  of  v/YLWL  allowed 

depending  upon  the  value  of  C  .   These  limits  are  presented  in  Table  V. 

b 

TABLE  V 
Limits  of  V/YL'/L  Alloyed  "by  Program 


Range  of  C 
b 

(v/Vlwl)   . 
v   '           'nun 

(v/Vlwl) 

max 

.600  <  c    <  .612 

"b  " 

•^ 

1.10 

.612  <  cb  <  .625 

A5 

.90 

.625  <  ao  <  .-J25 

Ao 

.90 

.725  <  cb  <  .775 

•  35 

.675 

•775  <  Cb  £  .000 

•  35 

.8co 

.26- 


Within  the  ranges  of  independent  variables  stipulated  above,  one  can 
expect  values  of  C  and  PC  from  equations  (7)  and  (8)  to  agree  with  the 
published  results  of  Series  60  to  within  Yp,   with  consistent  reliability  . 
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IV.  DESCRIPTION  OF  TEE  SEP-5HP  COMPUTER  PROGRAM  FOR  SERIES  60 

The  program  for  calculating  EHP  and  SH?  using  the  Series  60  model 

test  results  vas  written  in  the  FORTRAN  IV  language.  A  flow  chart  for 

the  program,  a  listing  of  the  FORTRAN  statements,  and  the  900  elements 

of  the  matrices  [A   I  and  |3  , 1  mentioned  in  Chanter  II  are  included  in 
L  nlJ     l  nlJ 

Appendix  A.  The  following  is  a  general  description  of  the  methods  of 
calculation  used.  Detailed  instructions  for  the  use  of  the  program  are 
contained  in  Appendix  C. 

1.  Inputs 

First  the  elements  of  matrices  (a  A  and  |3  ,1  are  read  in  as  data. 

1  nl J     l  iuJ 

The  number  of  ships  to  be  calculated  (NJOBS)  and  an  identifying  alphameric 
statement  is  then  read.   It  takes  three  cards  to  describe  each  ship  to 
be  calculated.   The  first  card  describes  the  input  mode  (lIIOPT)  for  the 
ship  to  be  calculated.   The  operator  is  afforded  an  option  of  one  of  two 
possible  input  modes  for  each  ship. 

a.  Mode  I 

Entry  is  made  with  L,  B,  H,  C  .  V^.,  Vmin,  Vincrement,  ACf ,  Position 
of  LC3  from  £§  ,  S,  *f ,   and  kinematic  viscosity. 
b.  Mode  II 

Entry  is  made  with  L,  B/H,  A,  Cp,  V^,  Vmin,  Vincrement,  ACf,  Position 
of  LCB  from  ^  ,  S,  ~£>,   and  kinematic  viscosity. 

2.  Calculation  of  needed  parameters 

Using  the  input  information,  the  program  then  calculates  the  remaining 
parameters  required  for  executing  the  remainder  of  the  program.   For  example, 
if  the  entry  is  made  using  Mode  I,  the  program  calculates  £,  3/h,  C  using 
equation  (l),  the  Series  oO  equivalent  LWL  using  equation  (2),  the  volume  of 
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displacement  ^  ,  midship  coefficient  C  ,  and  L/B. 

3.   Normalizing  hull  proportion  parameters 

The  program  calculates  the  normalized  hull  proportion  parameters 

for  L/B,  B/H,  and  C  using  the  eauations  in  Table  II. 

b 


k.      Trap  to  make  sure  V/VU-rL   is  within  range  of  model  test 

The  program  calculates  v/VTJ/L  maximum  and  minimum  and  ensures  that 
these  values  are  -within  the  range  of  validity  given  in  Table  V.   If  they 
are  not,  the  program  alters  V  ..  and  V    accordingly,  until  they  are 
within  the  range. 

5.  Calculation  of  number  of  increments  in  V  to  be  calculated 

Using  VKiax,    v  .  ,  and  V^ncrerrien-t,  the  number  of  different  speeds  (m) 
at  which  EHP  and  3H?  are  desired  is  calculated.  There  is  a  maximum  of 
100  different  speeds  allowed  for  each  ship.   If  this  number  exceeds  100, 
the  program  keeps  doubling  Vincremen-{;   until  "m"  is  100  or  less.   If  the 
operator  desires  EHP-SKP  information  for  only  one  speed,  m  is  set  at  one. 

6.  Calculation  of  V's,  v/YLWL's  and  normalizing  V/YlWL 

The  V's  and  v/Vlv/L's  are  calculated  based  on  the  VU]±n      ,  Vincrenient 
LWL,  and  m  calculated  above.   The  v/VlUL's  are  then  normalized  using  the 
equation  in  Table  II.    If  V^^v  =  V^q  and  v"increment  =  0.,  only  one  speed 
is  used. 

7.  .  Calculation  of  C^  for  a  UOO  foot  ship  and  PC  for  a  600  foot  ship 

The  C  for  a  ship  of  L  =  ^00  ft.  (C,  )  and  PC  for  a  ship  of  L  =  600 
ft.  are  calculated  for  each  v/VLWL  -using  the  normalized  parameters  determined 
above,  and  equations  (7)  and  (8)  respectively.   Subroutine  CEF  is  used  to 
*  See.  program  listing  in  Appendix  A. 
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evaluate  the  Chebyshev  Polynonials  for  these  equations. 


8 .   Calculation  of  C^  for  a  400  foot  ship  and  C*.   for  a  ship  of  length  L 


The  Reynolds  Numbers  (Re)  for  the  L  =  ^00  ft.  ship  (Re„J  and  for  the 

M 

ship  of  length  L  (Re  )  are  calculated  for  each  v/YlWL  by  using  equations 
(11)  and  (12). 

*     406.75  x  vs 
ReM   =  : S  _  (11) 


M 


1.2817  x  10-5 


LWL  x  v0  (l2) 

ReS    =  5 

u 

The  C  's  for  the  ^00  ft.  ship  (C«  )  and  shin  of  length  L  (Cp„)  are 
f  fM  fS 

calculated  using  the  ATTC  19^-7  correlation  line  given  in  equation  (l3)« 
Cf"5  =  4.132  log10  (Re  x  Cf )  (13) 


In  the  calculation  of  the  C  's,  it  is  necessary  for  the  program  to  use  an 
iteration  process,  since  C  appears  on  both  sides  of  equation  (l3)«   The 
initial  estimate  for  this  process  is  made  by  using  the  equation  for  the 
International  Towing  Tank  Conference  (ITTC)  line,  equation  (l4). 


Q.     ,      0.075  (^) 

f      (log10  Re  -  2)d 
Equation  (13)  is  recycled  for  50  times,  or  until  the  C  values  converge  to 

vithin  a  tolerance  of  5.0  x  10  "  '  .  Although  this  should  not  occur,  if 

formula  (13)  does  not  converge,  the  program  prints  the  error  signal 

*  Note  that  LWL  =  ^Go.75  for  a  ship  of  400  ft.  length  between  perpendiculars, 
and  that  kinematic  viscosity  is  for  salt  water  at  59^* 
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"SCHOENHER  NO  CONVERGE"  with  the  speed  at  which  the  trouble  was  experienced. 

9.  Calculation  of  wetted  surface  (S)  if  not  input 

Ideally,  the  wetted  surface  (S),  in  square  feet,  is  fed  into  the 
computer  as  input.   However,  if  the  wetted  surface  is  not  known  at  that 
time,  the  program  will  calculate  an  approximate  wetted  surface  using 
equation  (l5)« 

,  .329   ,  .203     -.0609     .666  ,       ' 

S  =  2.65*  (L/B)     (B/H)     (CD)       (V)  (15) 

This  equation  for  approximating  wetted  surface  for  the  Series  60  was 
developed  by  D.  Khoushy  in  reference  (2).   This  formula  was  tested  exten- 
sively and  conroarpd  to  the  contours  for  wetted  surface  calculation,  Figures 
B-12U  through  B-126,  in  reference  (l).   It  was  found  S..:.t   this  formula 
adequately  represents  the  wetted  surface  contours  of  the  Series  60,  with 
the  maximum  error  rarely  exceeding  1.2$.   Since  wetted  surface  is  a  factor 
in  calculating  resistance,  EHP  and  SEP,  if  wetted  surface  is  not  an  input 
parameter  to  the  program, Dne  can  expect  errors  in  EEP  and  SEP  due  to  the 
use  of  equation  (l5)>  in  addition  to  the  errors  to  be  expected  due  to  the 
use  of  equations  (7)  and  (8)  as  discussed  in  Chapter  III. 

10.  Calculation  of  C  for  the  ship  of  length  L,  R,  EEP,  and  SEP 


The  C  for  the  ship  of  length  L,  (C+  ),  is  calculated  using  equation 

t  b 


(16). 


C+   a  C+  -  CL  -  0.000^  +  C-  +  ac^  (16) 

tS      tM    fM  fS     f 


Then  the  remaining  calculations  are  trivial,  namely: 
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R  =  1/2  >f   3  v  2  C 


EHP  =  R  v 
s 


550 

SEP  as  EHP 
PC 

vhere  PC  is  for  a  ship  of  L  =  600  ft. 

11.  Calculation  of  LCB  position  from  ^  as  a  percentage  of  L 

The  LCB  position  from  te  as  a  percentage  of  length  between  perpen- 
diculars is  calculated  and  presented  in  the  output.   It  should  be  remembered 
here  that  corrections  to  R,  EHP,  and  SEP  for  variation  in  LCB  position  from 
those  of  the  parent  model  hull  forms  are  not  built  into  the  program.   Graphs 
for  assistance  in  making  these  corrections,  basea  on  Tables  ka,  through  53 
of  reference  (l),  are  included  in  Appendix  B.   The  LCB  position  from  $£ 
as  a  percentage  of  length  between  perpendiculars  is  calculated  for  entry 
to  these  graphs. 

12.  Output 

The  output  for  each  ship  calculated  includes  C  ,  L/B,  B/H,  C  ,  Cm, 

length  between  perpendiculars  (LBP)  in  feet,  LWL  in  feet,  3  in  feet,  H  in 

feet,  displacement  (DIS?)  in  tons,  volume  of  displacement  (VOL)  in  cubic 

feet,  vetted  surface  (WSUP.F)  in  square  feet,  AC-,  vater  density  ^f  in 

lb .  se c  ,  kinematic  viscosity  (GIU)  in  ft  ,  LCB  position  as  a  percentage 
ttJ  sec 

of  LBP,  and  a  tabulation  of  V,  v/YlTTl,  R,  EHP,  PC,  and  SEP  for  each  speed 

requested  in  the  input.  A  statement  is  printed  under  this  tabulation  vhich 
reminds  the  user  that  the  values  are  not  corrected  for  variation  in  LCB 
position  as  explained  above. 
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V.  EXAMPLES  IK  USE  OF  THE  EHP-SHP  COMPUTER  FROGPAH  FOR  KNOWN  SHIPS 

The  final  EHP-SHP  computer  .program  was  tested  with  the  parameters 

of  six  ships  whose  EHP-SHP  characteristics  are  known.  As  a  first  test, 

the  Series  60  equivalent  of  the  ship  SCHUYLER  OTIS  BLAND  was  chosen, 

since  this  was  the  ship  used  in  Appendix  D  of  reference  (l)  in  illustrating 

a  numerical  example  of  the  use  of  the  Series  60  charts.   It  was  also 

chosen  since  it  illustrates  the  ability  of  the  program  to  interpolate  for 

intermdeiate  values  of  C   (0.651)  and  l/b  (6.82)  and  to  extrapolate 

b 

outside  the  range  of  the  Series  60  in  B/H  (2.kkk). 

As  a  second  test,  the  five  parent  models  of  the  Series  60,  with 
L  =  600  feet,  were  chosen  since  complete  information  on  their  EHP  and  SHP 
"behavior  is  presented  5n  Tables  27  through  31  in  reference  (l),  and  they 
cover  a  variety  of  hull  proportion  combinations. 

The  results  of  these  tests  are  discussed  below. 

1.   Test  using  the  ship  SCHUYLER  OTIS  BLAND 

As  a  first  test  on  this  ship,  input  Mode  I  was  used  with  the  wetted 
surface  area  of  3999^*0  square  feet  as  input.  As  a  result,  the  computer 
output,  Table  VI,  for  EHP  was  a  measure  of  the  accuracy  of  equation  (7)« 
In  comparing  the  computer  program's  EHP  results  with  those  resulting  from 
the  hand  calculations  carried  out  in  Table  D-3  of  reference  (l),  it  was 
found  that  the  maximum  error  experienced  in  EHP  was  +  2.70$  at  v/YlWL  =  .85, 
while  the  average  error  experienced  was  +  1.19$  •  Although  the  average  error 
of  1.19$  is  fairly  gratifying,  the  maximum  error  certainly  is  not.  However, 
referring  to  Table  IV  again  we  see  that  in  using  the  hull  proportions  of  the 
SCHUYLER  OTIS  BLAND  we  place  ourselves  in  the  neighborhood  of  parent  model 
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U2l8  vhere  the  maximum  error  of  2.7^$  was  experienced  in  the  test  of 
equation  (7)« 

The  SCHUYLER  OTIS  BLAH)  was  also  tested  without  wetted  surface  as 
input.  The  results  of  this  run  are  presented  in  Table  VII.  We  can  see 
that  the  wetted  surface  calculated  by  equation  (15),  39651*5^  square  feet, 
is  in  error  by  -.858$.  The  error  in  the  calculation  of  wetted  surface  area 
tends  .to  reduce  the  error  in  EHP  in  this  instance. 

One  must  keep  in  mind  that  we  are  comparing  computer  output  with 
values  derived  from  hand  calculations,  in  which  design  charts  were  read  with 
an  unknown  accuracy. 
2.   Test  using  the  five  parent  models  of  Series  60  with  L  =  600  feet 

The  results  of  the  computer  run  with  the  five  parent  node Is  of 
Series  60  with  L  =  600  feet  are  presented  in  Tables  VIE  through  XII.   In 
these  cases,  the  wetted  surface  area  of  the  ships  was  calculated  by  equation 
(15),  since  it  was  not  known.  The  maximum  and  average  percent  errors  in 
EHP  and  3H?  experienced  in  this  test  are  presented  in  Table  XIII.  Although 
at  some  points  the  errors  exceeded  2$,  in  no  case  did  the  average  error 
exceed  1$.   These  errors  are  probably  due  mainly  to  the  inna curacy  of  the 
equation  used  for  wetted  surfa.ee  area. 
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VI.  CONCLUSIONS  ALT)  RECO^-SNDATIONS 

Based  upon  the  results  of  tests  run  to  date,  one  can  conclude  that 
the  EHP-SHP  computer  program  presented  herein  can  he  relied  upon  to 
estimate  the  EHP-SHP  requirements  of  single-screw  merchant  ships,  falling 
vithin  the  range  of  hull  proportions  of  the  Series  60,  to  a  tolerance  in  the 
neighborhood  of  1$  on  the  average,  and  with  great  certainty,  to  within  a 
tolerance  of  less  than  3$  f°r  any  given  speed. 

The  weak  link  in  the  EHP-SHP  program  appears  to  be  the  calculation 

of  vetted  surface  area.   It  is  recommended  that  the  representation  of  the 

Series  60  wetted  surface  contours  he  further  investigated,  to  determine 

whether  equation  (15)  can  he  modified  and  improved. 

It  is  further  recommended  that  the  elements  of  matrices  f  A  "]  and 

k  nl' 

[  'Bnj\   be  punched  in  the  form  of  FORTRAN  data  statements,  and  included  in  the 
body  of  the  main  program.   This  will  greatly  facilitate  its  conversion  to 
a  subroutine  for  inclusion  in  a  larger  program  if  it  were  so  desired. 

It  is  also  recommended  that  the  program  be  extensively  tested  with 
existing  ships,  to  further  determine  its  accuracy. 
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1.  Flov  chart  for  EHP-SHP  program.  hj 

2.  Listing  of  EHP-SHP  program  and  S/R  CEF.  ^9 

3.  Listing-  of  the  elements  of  matrices      |a    1    and  [b     |    .  53 


1,     Flow  Chart  for  EHP-SKP  Program. 


I       START        J 


READ: 


V>  Vx*  ?JOBS, 
ID  Information    / 


\ 


1 


(+■> 


0 


<-) 


1 


READ: 

ISP,  3,  H,   Cb 
VHIN,  VMAX,  VIKCi 


Calculate: 

A  *  3/H,  C, 


READ : 
L3P,  B/H,  A  ,  CD 
VMIN,  VMAX,  VINCy 
i 


READ : 
ACf,   LCB 
S,  RHO,  0N17 


Calculate: 


Calculate: 

EtfL,  \7  ,  C^,   L/3 

Normalize: 

1/3,  B/H,   Cb 


A 


)etermine: 
VLKAX,  VLMIT; 
from  C--, 


^ U!^X/VWL  :  VLM&xV 
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r~  C 
^7" 


Calculate: 

V  =7>I  II I  *  ?  x  VINC 

Nbrmalizpe  Vy  yUTl 


P=  P  +  1 


Call  S/fe  CEF 


x3   u, 


1 


Calculate: 

C^  ,   PC  for  each  V/"|/u7 


Cslcu3  ate: 

RfeK,  Re3,  CfM,  Cfs 


Calculate: 

CU_,  R,  EHP,  SHP 


% 


K3 


MJ03   :  NJOB3 
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2.      Listing  of  EHP-SHP  program  and  S/R  CEF 


1 
2 

3 
4 
5 

6 


PROGRAM    TOCALC    EHP,SHP    USING 
DI  4ENSION    A( 10,^5) ,B( 10,45) 
DIMENSION    VS("l6d)'f  VL(  106)"tV(lC 
DI  4ENS  [0  1    PSII  L'1  )_,CHIUC)  ,VPHI 
DIMENSION    CFM{  100  )  ,"cf~S  (  100  J  ,R( 
FORMAT     (5F16.71 

(  14*  1  -)  A  +  ) 

(  7  F 1  n  .  2  ) 

( 4F10.6, £10.6) 


SERIES  30  MODEL  TcST  RESULTS 


G  )  ,  A  V  L  (  1 0 
(  1 0 ) , C  TM ( 

100) ,  EHP{ 


),3VL(1C),PHIC1Q) 

100)  ,  PCM  100)  ,  CFS(  13C 
100)  ,SHP(  100) ,QQ(20) 


FORMA r 
F  0  R  M  \  T 
FUR  4AT 
FOR  4  AT 
FOR  '  \T 


( 


SCHOENHER 


NO  CONV  E 
2  X ,  5  m  l 


l4H0Cb  =  t  F  6 • 3 » 
IF 6.  3t  2X_,4H  CM=»_F6.3J 

7  FORMAT  (5H  Lyp=,F3.2,3H  FT,2X, 
iFd^2,3H  FT_,2X,3H  H=,F8.2,3H  FT 

8  FORMAT  (  6H  0  I SP= , F 10 . 2 , 5H  TON 
17H  rtSJRF=,_F10.2_,oH  FT**2) 

9  FORMAT  (1H0,2X,7:W  (  KTS  )',  3X~,  7H 
14X,2HPC,7X, 3HSHP  (HP) J 

10  FORMA*  (~F9.~3VF10.4»Fi4.2~,F13~.2 
500  FORMAT  (_fH  Lu3  IS,  F6_.  3,  33. 1  PER 
50  1  FORMAT  (21H  LOO  IS  AT  AMID SHIP 

502  F0R_MAT_j(_7_H_.J._C3  IS,C6.3,33H  PER 

503  FORMAT  ( 14) 

534  FOR  4AJ  (  1H1  ,  14,  L9A4) 

505  FORMAT  (  1GH  DELTA  CF=,F1C.7,2X 
i_2X,  5H  GNU=,E1_6.9_,1CH  FT**2/SE 

506  ^Oa-IAT  ("74H0"  ***  ABOVE  VALUES 
IOF  LO  j  POSITION  ***) 


RG£_.   V  S=  t  rb  .  3  ) 

/fcJ=,F6.3,2X»5H  B7h=»F6.3,2X,4H 


5H  LWL=,F 

) 

S  ,  2  X  ,  5  H  V 


B.2,3H  F  r ,2X»3M  B=i 
0L=,F12.2,6H  FT**3'f2X 
■J/'   TL  -.LroX,  7HR  (  LBS  )  ,  5X,  6HEHP 


(HP) 


,  F  8 . 4 ,  F  1'3 
CENT  Oc  L 
S.  ) 
C  h  N  T  0  F  L 


,  5  H  R  H  0  =  i 

AR'E   NOT   C 


.2) 

FjP  AFT  OF  AMIDSHIPS.  ) 

3P  FWD  OF  AM  IDSHIPS.  ) 

F10.7,13H  LB*SlC**2/lF 

O'RRECTcD  FOR  ANY  VA<I 


T  *  *  ^ 

at  i : 


200 


201 


9rC 


READ  COEFF  OF  POL  FOR  CTM  AtlJ    PCM 
DO  2rO  0=1, 10 

RcAj     (3,1)  (aTn,L  )  ,L=1,45) 
CONTINUE 
00  2C  1  i=l,  1 

READ  (5,1)  (0 (  i,L  )  ,  L=l,45) 
CON I  I OOc 

READ  NUMBER  OF  SHIPS  AND  10  INFO 
READ  (3,2)  NJU3S,  (00(L)  ,  L=l,  i  i] 
OU    1000    NJOB=l,NJOBS 

READ    INPUf~AN0    CALC    NEEDED    PARAMETER'S    FOR    SHI=> 
REAO    ( 5,503 )     I  iOPT 
IF     (  IruPT  J     900, 900, hCl 

s.:AO(t,3)     XL6P,     3H,     DEL,CP,  VMIN,  VM_AX,  tfINC 
C  f=  .  '•  357*(  OP  )**"2    +     7925*CP 
H=SQRT  (  (  35_.*DEL  )  /  (C3*  <LjJP*BH)  ) 
B*M=BH*H 
GU     TO    903 


901     READ  (5,3)     aL.P,     3M",     H,    C 
DEL=C5*XL  3P*BM*H/35. 


VMIN,     VfrUX 


-1,9- 


(03 


600 


ell 


802 


8C  i- 
II 

12 

13 
14 

15 


16 


CP=(-.925    +    SQRTC.  925**2    +    4.*  .  0375*CB  )  )  /  (  2.  *  .  0375  ) 

READ    l.5f-4)pCjF»XLC.-ii  WSUitF»RH0»G»\U  ._     __.__ 

XLWL=40&.75*XLQP/4CC 

^/OL-B5-*DEL 

CM=C3/CP 

NOR  -IAL  I  Zf     JA;U  -IE  Ft  *S 

U=i^ 7  . 0  -  *  -B L  J  / 1  .  5 

Z=-6.C     *    2.C*BH 

X  =  -7. 0    *■    lQ.0*Cd— 

TRAP    TO     -TAKE    SjxE    VL     IS    WITHIN    RANGE     Jc    MODEL    TEST    RESULTS 

IF     LC3.GE.C..6G.ANQ.C3.LE.G.612)- GO    TO    830      - -     - 

IF     (C3.GT.0 .612.AND.C6.LE.0.625 )    GO    TO    801 
IF     CC3.GI.C.625.AND.C8.LE. 0.725 1    00    TO    802- 
IF     (C3.GT.0.725.ANO.C3.LE.0.775)    Go     TO    303 
IF     (.Ca»GT--0.775«ANO»C-B.LE. 0.300)    GO     TO    804 
VLO  [N=.45 
VL.HAX=1.-1 

go   ru   n 

VIM IH- .45 

V  L "'  \  X  =  .  ^  Q 

Gu-ian  .____._. 


tfLHAXs .  r 

^■J    TO     1  1 

VL..r:=.  3  5  — _ 

■J  L  '  A  X  =  .  j  7  5 

5U-IU11  

7L  IIN=.  35 
\/-L  M  A  a.  =  *.  3  0 
I  F     (  V  •'  i    1/  (  XL  *L**  .  5  )  - V  L ' !  I  N  )  1  2 ,  1 3  ,  1  3 


LrtIfi*lXLWL**»5] 


v.-u 

[F     (VLiMAX-VHAX/(XLWL**.5n  14,  15,  15 

■/  '1  iX=VLMAX*  l  X.L*IL**»5  1 

CALC  OF  -i.  IF  VINC=0.,M=1  THROUGHOUT  PROGRAM 

IF    IV  I  iC.'jI .  :  .  )    Go    TO     L5 

M  =  l 

iMIN>V  4AX  .    . 

00    TO    20 

y/M=  CV.'4AX-VMI:\I1/V1MC    +    1.5 

■1  =  7  0 

Et=a    +    1  ______ 

10\E    SURE      ;      iOT    GREATER    THAN     LOG 

IF     Lrt.LS~.iCQ)     OJ     fj    20 

7  I  4C  =  2.  ;*v  J   iC 

GO    TO     Id 

_ALC     JF     SHIPS    v,     v     DVcK    Of    L  ■.  L ,      IURM-.LIZE     v    UVER    ST    LWL 


-50- 


-2t- 


1-2 


32 


36 

40 


41 
4.5 

50 


J  J '  3  J  [ 
V S I  I)  =  v'' 

V  L  (  I  )  =  V 

V  <  I  )  =  2  , 
P=P+I.O 
Cil^TINU 
C  A  L  C  J  F 
CALL  CE 
CALL  ce 
CALL  CE 
DO  3  2  N 
AvL  CM)  = 
3  VL  (  N )  = 
CQNT LMU 
DO  35  K 

.DO  35 
DC  3d 
L=J  + 
DO  35 
AVL  1 14  j  = 
3VL(.M)  = 

DO  40  I 
CALL  CE 
G  TM ( I M  ) 
BCM-l  I  U 

DO  35  N 

cr.K  in 

PCH(  EM  J 

CT.M  I  -1) 
CON  T I N U 
C  AL  C  iiF 
00  50  I 
.REt4.-l.V-S 
R  E  S  =  (  V  S 
Cr,;(  I  )  = 
CFS( I j = 
DO  4d  K. 
DO  I1Y  =  C 
3QGUS=C 
C  F  H  I  )  = 
C  F  S 1 1 )  = 

Ir( ABSt 

1ELA3S( 

CO  ^T  loo 
WRITE  ( 
CO  iTINU 


=   L  1    " 
■1  I     J       * 

s  ( r )  / 

5*VL  ( 


CT'  S 

F  (  5  ,  X 
F  t  3 ,  U 
F13,Z 

=  L ,  1  C 

-\        O 

—    •  \_* 

).  ] 

c 

=  1,5 

=  1,3 
=  1,3 
*  (  I  -  I 
=  1,10 
aVL  t  N 
riVL  IN 

-1=1,  M 
r  {   I  C  , 

=0  .  c 

=.  (j  «-0 
=  1,10 
=  CTMl 
=  PCM( 
=CTH.( 
E 

CFIS 
=  1 ,  M 
111*1 
(  I  )  *  1 
0  .  0  7  5 
:■ .'  0  7  5 

=  1,5: 

F  -!  (  I  ) 
F.S(.L) 

(2.30 
(  2.30 

JJ.-i.-lY 


P*V-INC 
(  XL  *L**.5.) 

I  )  -  I .  d  7 5 


A  10    JC  '  S    FOR 
,  C  .4 1.  J 
,PHIJ 
,PSIJ 


4  00    FT    MODEL 


)     *     3*(K-1) 

>+A{  -. .  Li*?HLU  1*?SI  1 1 )  *CH  I  t  K.) 

)+3(n,d*phi  { j)  --  „:;;,;  c:m  * » 


/(  1-1)  r-VPHll 


LM)  +-AVL  INI*VPHI  U\\. 
IH)+BVL(N)*VPHI (  i) 
I  0)-:VL  (  UU*VL  (  LMMeG...QG.l 


|-J.<  400 


FT  MODEL  AND  UJ.^  OH  IP  LENGTH 


1  /c-05) 


.6839*406.  75J/CL. 2 
,6339*XLWL ) /( GNJi 
/  (  (  (  ALOC-tRE  »U/2  .30  25351 1-2.0) 
/U  (  ALOG(RES)/2.302  5ti51)-2.0) 


**2 

a,  »•-   "\ 


2) 
2) 


25  851/(4. 132* A LOG (REM*CFM( I ) ) 
25.351/14.  132*AL  J0(  Kc.S*CFS(  I  )  ) 
-CFM(  I  )  )-5.  ">E-C  7)  41, 4 1,45 
50GUS-CF  S I  I  )  )  -  5  . '  t - G  Z )  5 0  ,  50.,  45 


)  )* 
)  )  * 


-:2 
*2 


b,5) 
E 


VS(  I  ) 
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3~: 


.__3G1 


.30 
8 


CAL 
IF 

^  WSJ 

CAL 

CTS 
R4-I 

SHP 

>0  CON 
ARl 
rtRI 
elRI 

ARl 

MRl 

CAL 
PLC 
IF( 

Gu 
aiRI 

GO 

.-.  R  I 
rt  R  I 
-  OU 
rfRI 


b  Ur 
(  rtSU 
RF  =  2 
C  OF 
60  I 
(  I)  = 
L=0- 
(  i  )  = 
I  [  )  = 
T  [NU 
Tc  ( 
TE 
TE 
TE 
Tc 
C 

J  =  A  Ji 
X  L  C  B 
TE     ( 


t>n 


TE 
TJ 
Tc 
TE 
70 
TE 


LQO! 


wRI 

CON 

s  ro 

ENQ 


TI..J 
TE  ( 
IINU 
p 


0 

=  1 
CT 
5* 
(R 

EH 
E 

a» 
6  ♦ 
6  » 
6  , 
o  » 
SI 
S( 
)3 
fa  > 

6  , 
0 

6  » 
fa  » 

=  1 
a , 
E 

fa  t 
E 


SURF     E-F    NOT    GIVE  n 

.GT  .0. J     GO    TO    5  2 

54*ldL**.329)  ndH**. 206 ) *l  U/CS**»06091* UOL**. 666) 

j  T  P  U  T 

(  I ) -CFM(  I  )-0.0004+CFS(  I  )+00F 
EtHQ+WSURF:*  (  (  y/S(  I)  *  l.ad*9  )  **2  )  -CTS  (  I  ) 
(I)*VS( I )*1.6889)/550. 
?(  I  )/PCM(  I  ) 

5(K.)     NJ06S,     LQQIL1,L-1*193  .       .  ._ 


) 


o  i  8  L  ,   i  H 


»  w  r  t  o . 


7)  XL3P.».XLWLt3rttH 

8)  DEL,  tfOL,  tfSURF 
505  J  DCF  » -;hU  »  GNU 

T  OF  LC->  AS  PERCENT  Lb 
XLCH*L0Q./XL8PJ 


500  J 

50  11 


Lca 


9) 

10)  VS(  I  )  , VL (  I  )  ,R(  I  )  , EHP(  [J ,PCM( I) ,  SHP( I  ) 

506) 


i-J5R-)JTI  NE  CEr(  1-^,0,  CEi) 

01  PENSION  CE7 (  10) 
-  C  E  V  I  1  )  =  1  . 

CEV< 2)=J 
—00    l_I=3vIQ 
1       CEV  (  I  )  =2  .*U*CE V  (  1-1 ) -CE\/U -2) 

R.ETOilN 

END 
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APPENDIX  B 


GRAPHS  FOR  USE  IK  CORRECTION 
OF  EHP  AMD  SHP  FOR  VARIATION  IN  LCB  POSITION 
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IMSTRUCIICTJ  E3  USE   0?  THE 
EHP-3HP  PROGRAM  FOR   SERIES  60 
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1.  Input  deck  for  matrices  fA  A      and  [3  -,*] 

The  input  of  matrices  f.A.  ,^  and  fs  "1  ,  a  total  of  900  elements, 
is  accomplished  by  reading  the:.!  from  a  deck  of  lSo  cards,  with  5 
elements  on  each  card,  arranged  in  Format  5El6»7«   See  Appendix  A  for 
their  proper  arrangement. 

Total  no.  of  cards  =  180 

2.  Identification  card 

In  columns  1  -  k,    (ITJ0B3),   Punch  the  integer  number  of  ships  to  be 

calculated  in  this  particular  run,  with  the 
units  digit  of  the  integer  in  column  h. 
Format  1^. 

In  columns  5  -  80,  (QQ)     Punch  any  identification  information  desired. 

Example  :   "SCHUYLER  OTIS  B1AIJD   TEST  RUN.  " 
Format  19A^. 

Total  no.  of  cards  =  1 

3*   Input  for  Sh'.ps 

a.   If  input  Mode* I  is  being  used  for  a  ship,  punch  the  following  3.  cards 

for  each  ship: 

Card  1 

In  columns  1  -  k,   .(iNOPT),  Punch  a  positive  integer.  Example:  "+5". 

Format  I^v. 

Card  2 

In  columns  1  -  10,  (XLBP),  Punch  the  length  between  perpendiculars  in  feet. 

Example:   "  "4  00.0"  Format  Fib. 2. 

In  columns  11  -  20,  (Bl-l),   Punch  the  maximum  beam  at  the  load  water  line  in 

feet.   Example:  "66.0".   Format  F10. 2. 

In  columns  21  -  30,  (H),   Punch  the  draft  in  feet.  Example:  "27-0". 

Format-  F10.2. 

In  columns  31  -   ho,    (C3),      Punch  the  block  coefficient.      Example:    "O.65I". 

Format  F10.2. 

In  columns  kl  -   50,  (vHErl),  Punch  the  lowest  jpeed  in  knots  for  which  you 

desire  information  on  EHP  and  3HP.  Example : 
"3.557".  Format  F10.2 
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In  columns  51  -  60,  (Vi-AX),   Punch  the  highest  speed  in  knots  for  which 

you  desire  information  on  EKP  and  SHP. 
Example:   "19.25".  Format  F10.2. 

In  columns  6l  -  70,  (VIHC)j   Punch  the  increment  of  speed  at  which  you  wish 

information  on  EHP  and  SHP.  Example:  "1.07" 

Format  F10.2. 


Card  3 

In  columns  1  -  10,  (DCF), 


Punch  the  frictional  ship  correlation  factor 
you  desire  to  use.   Example:  ".OOOU". 
Format  F10.2. 


In  columns  11  -  20,  (XLCB),   Punch  the  position  of  LCB  from  ^  in  feet, 

with  a  +  sign  if  forward  and  a  -  sign  if 
aft  of  ££  .   Example:   "-  k.O".   Format  F10.2. 

In  columns  21  -  30,  (WSURF),  Punch  the  wetted  surface  area  in  square  feet. 

Example:   "3999^.0''.  Format  F10.2.   If  wetted 
surface  is  not  known,  punch  a  "0.0"  in  this 
space  and  the  program  will  calculate  it. 


In  columns  51  -  kO,    (REO), 


In  columns  kl  -   50,  (GNU), 


2/   * 

Punch  the  density  of  water  in  lb.  sec  /ft. 

at  the  temperature  desired.  Example:  "1.9905"' 
Format  F10.2. 

Punch  the  kinematic. viscosity  in  ft.  /sec. 
Example:  "1.2817E-05",  for  1.2817  x  10"5. 
Format  E10.2. 


Note :   If  information  on  EHP  and  SHP  is  desired  for  only  one  speed  for  the 
ship,  e.g.  service  speed,  on  card  2  punch  this  speed  in  the  spaces  for 
VMIN  and  VMAX,  and  punch  "0.0"  in  the  space  for  VINC. 

Total  no.  of  cards  =  NJ0BS  x  3. 

b.   If  input  Mode  II  is  being  used  for  a  ship,  follow  the  instructions 

for  Made  I  above  with  the  following  exceptions: 

Card  1 


In  columns  1  -  k,    (IN0PT), 


Punch  a  negative  integer.   Example:  "-5". 
Format  1^. 


Card  2 


In  columns  11  -  20,  (3H), 


Punch  the  beam  to  draft  ratio.  Example 
"2.5".  Format  F10.2 


In  columns  21  -  30,  (DSL),   Punch  the  displacement  in  tons.   Example: 

"10000.0".   Format  F10.2. 


-76- 


In  columns  31  -  ^0>    (CP),    Punch  the  prisiratic  coefficient.   Example 

".6l2".   Format  F10.2. 

Card  3 

Identical  to  Mode  I. 

Mode  I  and  Mode  II  input  ships  can  be  mixed  at  will  as  long  as  the 
cards  are  punched  as  above. 

Graphs  for  making  corrections  to  the  computer  output  for  variation 
in  LC3  position  are  contained  in  Appendix  B. 

k.      Loading  of  cards 

The  cards  are  loaded  in  the  following  order: 

1.  Main  EKP-SHP  Program  Deck 

2.  Subroutine  CEF 

3«      Input  deck  for  matrices     [A    "^  ancrB^j 
k.      The  identification  card 
5.      Input  cards  for  ships. 
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